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Abstract. With the first light of COROT, the preparation of KEPLER 
and the future helioseismology spatial projects such as GOLF-NG, a 
coherent picture of the evolution of rotating stars from their birth to 
their death is needed. We describe here the modelling of the macro- 
scopic transport of angular momentum and matter in stellar interiors 
that we have undertaken to reach this goal. First, we recall in detail 
the dynamical processes that are driving these mechanisms in rotating 
stars and the theoretical advances we have achieved. Then, we present 
our new results of numerical simulations which allow us to follow in 2D 
the secular hydrodynamics of rotating stars, assuming that anisotropic 
turbulence enforces a shellular rotation law. Finally, we show how this 
work is leading to a dynamical vision of the Hertzsprung-Russel dia- 
gram with the support of asteroseismology and helioseismology, seismic 
observables giving constraints on the modelling of the internal trans- 
port and mixing processes. In conclusion, we present the different 
processes that should be studied in the next future to improve our 
description of stellar radiation zones. 



1 CEA/DSM/DAPNIA/Service dAstrophysique, CEA/Saclay; AIM-Unite Mixte de Recherche 
CEA-CNRS-Universite Paris VII, UMR 7158; F-91191 Gif-sur-Yvette Cedex, France 

2 LUTH, Observatoire de Paris; 5 place Jules Janssen, F-92 195 Meudon Cedex, France 

3 Observatoire de Geneve; 51 chemin des Maillettes, CH-1290 Sauverny, Switzerland 

4 Institut dAstrophysique et de Geophysique, Universite de Liege; Allee du 6 Aout, 17 1/14 
(B5C), B-4000 Liege 1 (Sart Tilman), Belgium 

5 GRAAL, Universite de Montpellier II; Place Eugene Bataillon, F-34095 Montpellier Cedex 
05, France 

6 Institut d'Astronomie et dAstrophysique, Universite Libre de Bruxelles; Campus de la 
Plaine, Boulevard du Triomphe, CP 226, B-1050 Bruxelles, Belgium 

7 LATT, Observatoire Midi-Pyrenees; 14 Avenue Edouard Belin, F-31 400 Toulouse, France 



© EDP Sciences 2008 
DOI: (will be inserted later) 



2 



Title : will be set by the publisher 



1 Dynamics of stellar radiation zones and differential rotation 

Rotation, and more precisely the differential rotation, has a major impact on the 
internal dynamics of stars. 

First, as it is known from the theory of rotating stars, rotation induces some-large 
scale circulations, both in radiation and convection zones, which act to transport 
simultaneously the angular momentum, the chemicals but also the magnetic field 
by advection. In radiation zones, the large-scale circulation, which is called the 
meridional circulation, is due to the differential rotation, to the transport of angu- 
lar momentum and to the action of the perturbing forces, namely the centrifugal 
force and the Lorentz force (cf. Busse 1982, Zahn 1992, Maeder & Zahn 1998, Ga- 
raud 2002, Rieutord 2006). Next, the differential rotation induces hydrodynamical 
turbulence in radiative regions through various instabilities: the secular and the 
dynamical shear instabilities, the baroclinic and the multidiffusive instabilities. In 
the same way that the atmospheric turbulence in the terrestrial atmosphere, this 
hydrodynamical turbulence acts to reduce the gradients of angular velocity and of 
chemical composition; it is why it is modelled like a diffusive process (cf. Talon & 
Zahn 1997, Garaud 2001, Maeder 2003). On the other hand, rotation has a strong 
impact on the stellar magnetism. For example, it interacts with the turbulent 
convection in convective envelopes of solar-type stars (cf. Brun et al. 2004) to 
lead to a dynamo mechanism and, as it is expected from observations, to a cyclic 
magnetism. In radiation regions, it interacts with fossil magnetic fields where the 
secular torque of the Lorentz force and the magnetohydrodynamical instabilities 
such as the Tayler-Spruit instability and the multidiffusive magnetic instabilities 
have a strong impact on the transport of angular momentum and of chemicals (cf. 
Charbonneau & Mac Gregor 1993, Garaud 2002, Spruit 1999, Spruit 2002, Menou 
et al. 2004, Maeder & Meynet 2004, Eggenberger et al. 2005, Braithwaite & 
Spruit 2005, Braithwaite 2006, Brun & Zahn 2006). Finally, waves constitute the 
last transport process in single stars and they are also interacting with rotation. 
Internal waves, which are excited at the borders with convective zones, propagate 
inside radiation zones where they extract or deposit angular momentum where 
they are damped leading to a modification of the angular velocity profile and thus 
of the chemicals distribution (cf. Goldreich & Nicholson 1989, Talon et al. 2002, 
Talon & Charbonnel 2003-2004-2005, Rogers & Glatzmaier 2005). Note also that 
rotation modifies stellar winds and mass losses (cf. Maeder 1999) 
In close binary systems, where the companion could be a star as well as a planet, 
there are transfers of angular momentum between the star, its companion and the 
orbit due to the dissipation acting on flows induced by the tidal potential; that 
could be the equilibrium tide (cf. Zahn 1966) due to the hydrostatic adjustement 
of the star or the dynamical tide which is due to the tidal excitation of internal 
waves (cf. Zahn 1975). This dynamical evolution modifies the internal rotation of 
each component that have consequences on the properties of their internal trans- 
port. 

To conclude, all the processes, with which rotation interacts, transport angular 
momentum and matter that modify the internal angular velocity, the chemical 
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composition and the nucleosynthesis. Therefore, differential rotation has impera- 
tively to be taken into account to get a coherent picture of the internal dynamics 
and of the evolution of the stars. 

2 Advances in secular magnetohydrodynamics of stellar radiation zones 

First, we have studied the rotational transport where the angular momentum and 
the chemicals are transported by the meridional circulation and by the hydro- 
dynamical turbulence due to shear instabilities, namely the rotational transport 
of type I. We generalize its present modelling to treat simultaneously the bulk 
of radiation zones and their interfaces with convective zones, the tachoclines (cf. 
Mathis & Zahn 2004). Next, we derive a new prescription for the horizontal tur- 
bulent transport which is obtained from Couettc- Taylor laboratory experiments 
that allow to study turbulence in differentially rotating flows (cf. Mathis et al. 

2004) . However, the introduction of these two hydrodynamical mechanisms in 
stellar models fail to reproduce the observations of solar-type stars, because these 
have been slowed down by the wind during their evolution and hence the rotational 
processes are less efficient. Therefore, we consider the rotational transport of type 
II where the chemicals are still transported by the meridional circulation and the 
turbulence, but where the angular momentum is carried by another process, the 
two candidates being the magnetic field and the internal waves. So, we introduce 
the effect of a fossil magnetic field in a consistent way in taking into account the 
action of turbulence, differential rotation and meridional circulation on the field, 
but also its feed-back on momentum and heat transports (cf. Mathis & Zahn 

2005) . Next, we introduce in the modelling of the internal waves the effect of the 
Coriolis force (cf. Mathis & Zahn 2005) that allows to include the gravito-inertial 
waves in the description of the angular momentum transport. Finally, a coherent 
treatment of tidal processes has been derived. A detailed review of previous results 
obtained both on the rotational transport of type I and of type II is done by J. -P. 
Zahn in this volume. 

3 Modelling 

To get a coherent dynamical description of stellar radiation zones, the complete 
equations of magnetohydrodynamics have to be solved. To achieve this goal, as- 
sumptions are made. 



3.1 Main assumptions 

Here, secular magnetohydrodynamics and its consequences on stellar evolution 
are studied. Thus, secular time-scales associated to the nuclear evolution of stars 
are chosen. Moreover, low angular resolution (with an expansion in few spherical 
harmonics) is considered due to the turbulent transport behaviour in radiation 
zones (cf. next paragraph). Physical processes which have dynamical time-scales 
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and need a high angular resolution description, such as hydro or magnctohydro- 
dynamical instabilities and turbulence, are treated using prescriptions. This is the 
first step to achieve the highly multi-scales problem in time and in space of the 
dynamical stellar evolution that could not be yet studied with Direct Numerical 
Simulations. 



On the other hand, stellar radiation zones are stably stratified regions. Thus, 
the buoyancy force, which is the restoring force, acts to inhibit turbulent motions 
in the vertical direction. This leads to a strongly anisotropic turbulent transport 
for which that in the horizontal direction (on an isobar) is more efficient than that 
in the vertical one. Therefore, horizontal eddy-transport coefficients are larger 
than those in the vertical direction and the horizontal gradients of scalar fields 
such as rotation, temperature and chemical concentration are smaller than their 
vertical gradients. An horizontal expansion in few spherical harmonics is thus 
allowed, and we get for a scalar field, X: 

X(r,6,<p)=X(r,t) + 8X(r,6,<p,t) (3.1) 

where SX (r, 9, <p, t) = £i>o EL-i (r, t) Y t m (9, tp) with X (r, t) » X l m (r, t), 
X and SX being respectively its horizontal average on an isobar and its fluctua- 
tion, r, 6, ip are the classical spherical coordinates while t is the classical time. 



The dynamical equations such as the induction equation for the magnetic field 
or the Navier-Stockes equation are three-dimensional vectorial equations. Here, 
the aim is to study the secular magnetohydrodynamics of rotating stars and its 
consequences on stellar structure and evolution using stellar evolution codes which 
are mostly unidimcnsional. To couple them with transport equations, we expand 
vector fields such as macroscopic velocities or magnetic field in vectorial spherical 
harmonics (see Rieutord 1 987) like in stellar oscillations theory: 



OO I 

u{r, 6,<p,t) = J2Yl {< ^ *) [ y ™ V) + v l m (r, t) [v H >T {9, ip) 



1=0 m=—l 



-w l m (r,t) [v H A(y ; m (^)e r )]} 



where Vh — egde + e„ 



A- (3-2) 



These expansions allow to separate variables in transport equations. Thus, 
modal equations in r and t only are obtained and can be implemented directly in 
stellar evolution codes. 



3.2 Preliminary definitions 

The macroscopic velocity field is expanded as: 

V = r sin Oft (r, 9) e v + re r + U M (r, 9) + u (r, 0, ip, t) . (3.3) 

The first term, where fl (r, 9) is the internal angular velocity and e v is the unit 
vector in the <p direction, is the azimuthal velocity field associated to the differential 
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rotation. The second term corresponds to the radial lagrangian velocity due to 
the contractions and the dilatations of the star during its evolution, e r being the 
radial unit vector. The third term Um (r, 9) is the meridional circulation velocity 
field which has been presented before. Following the general method concerning 
the expansion of vector fields, it is projected on vectorial spherical harmonics: 

Um = {ui ( r ) p i (cos0)e r + Vi (r) dP ^ c » sg ) gA . (3 . 4) 

The anelastic approximation is adopted, thus filtering-out sonic waves, that is 
justified for the slow meridional circulation. Therefore, the continuity equation 

d t p + V ■ (pV^j = 0, where p is the density, becomes V • ^/oWm) = 0, that leads to 
the following relation between the othoradial functions, Vj, and the radial one, Up. 

( r ) — i(i+i) pr d ^ P d r U ^ ■ Finally, u(r,8,(p,t) is the velocity field of the internal 
waves. 



Next, the temperature, T, and the mean molecular weight, p, are respectively 
expanded as: 

T(r,9) = T(r)+5T(r,9) with 5T (r, 9) = Y [*; (r) T\ P t (cos 9) , (3.5) 

Z>2 



and p(r,6) = p (r) + Sp(r, 6) with dp (r, 9) = Y [A, (r) p] Pi (cos 6) ; (3.6) 

l>2 

T and ~p are their horizontal averages, ST and dp being their fluctuations and "f; 
and A; their relative fluctuations. 

Finally, the magnetic field is expanded using its divergence-free property: 

B (r, 9) = V A V A (£ P (r, 9) e r ) + V A (fr (r, 9) e r ) . (3.7) 

£p and £t are respectively the poloidal and the toroidal magnetic stream functions 
which are expanded in spherical harmonics: 

oo oo 

£p (r, 9) = Y, & M Y? (*) and (r, 9) = Y X < (r) r,° (0) . (3.8) 

z=i i=i 

Here, the mean axisymetric magnetic field is considered (m = 0). Therefore, the 
poloidal field, Bp (r,9), is in the meridional plane (e rj eg) while the toroidal one, 
Bt (r, 9), is purely azimuthal along e v . 

The different fields being now defined, we have to consider the transport equa- 
tions that have to be implemented in stellar evolution codes. 
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3.3 Transport equations system 

The first equation of transport we consider is the one related to the transport of 
the magnetic field, namely the equation of induction: 



3 ( B = VA (yAB) - VA (||t/||®VAb) 



(3.9) 



being the magnetic eddy-diffusivity tensor. Two advection-diffusion equations, 
for respectively £p and £t, are obtained from its spectral expansion in the vectorial 
spherical harmonics: 



dt 



rV A d:i (U U B) = Vh rA l 



and 



la 



Ila 



d r (VhdrXo) - Vvl (I + 1) 



(3.10) 



lb 



lib 



The terms la and lb, where Vaa-.i (Ui,E?j and 7Ad : ; Ui,B*j arc function of the 
differential rotation, the meridional circulation and the magnetic field, correspond 
to the advection of B by Um and to the production of its toroidal component 
through the shear of the differential rotation. Terms Ila and lib correspond to 
the turbulent ohmic diffusion, where the possibility of an anisotropic turbulent 
transport of the magnetic field is assumed, with r\ v and r\h being respectively the 
eddy-magnetic diffusivity in the vertical direction and in the horizontal one. The 
transport is studied from a lagrangian point of view due to the contractions and 
to the dilatations of the star during its evolution. The time-lagrangian derivative 
d/dt is defined by d/dt = dt + rd r . 

The second equation which has to be treated is the Navier-Stockes equation, 
in other words the equation of dynamics: 



d t V + 



(v-.v) 



V 



-VP- pV(j> + V- ||t|| + 



Mo v > 



KB. (3.11) 



P, <f> are respectively the pressure and the gravitational potential. ||r|| is the 
Reynolds stress tensor and the last term is the Lorentz force where p is the 
magnetic permeability of vaccum. Taking its azimuthal component and averaging 
over ip, the following advection-diffusion equation for the mean rotation rate on 
an isobar is obtained: 

P ii ^ ~h dr = h dr (^ 4 ^ n ) (?) h dr {Tj (r)] 



in 



IV 



(3.12) 
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Term I corresponds to the ad- 



where Sl(r) = (r, 6) sin 2 6d6/ f* sin 3 06.6. 
vection of angular momentum by the meridional circulation and term II is the 
diffusive term associated to the action of the shear-induced turbulence where v v 
is the eddy-viscosity in the vertical direction. These two first terms correspond to 
the rotational transport of type I. Term III is associated to the Lorentz torque, 
T^, and term IV corresponds to the transport by internal waves, Tj (r) being 
the associated mean angular momentum flux on an isobar. These two last terms 
correspond to the rotational transport of type II. The same type of equation is 
obtained for the differential rotation in latitude (cf. Mathis & Zahn 2004-2005). 

On the other hand, due to the long time-scales associated to the meridional 
circulation, dynamical terms in the meridional components of the Navier-Stockes 
equation are filtered, keeping only its hydrostatic terms. Taking the curl of the 
hydrostatic equation, we get the associated thermal-wind equation: 



<pAi-5*i = T -Vi (n,B 



(3.13) 



where ~§ is the horizontal average of gravity, the explicit form of T>i in function 
of £1 and B being given in Mathis & Zahn 2005. The more general equation of 
state is used (cf. Kippenhahn & Weigert 1990) where S = — (d\n p/ d\nT) p ^ and 
if = (d In p/d In p) p T . 

The last equation that has to be solved is the one for the macroscopic entropy, 
S, in other words the equation for the transport of heat: 



PT 



dtS + V-VS 



V- xVT +pe- V-F + J. 



(3.14) 



Its describes the transport of entropy by advection with taking into account the 
thermal diffusion (x is the thermal conductivity), the production of energy asso- 
ciated to the nuclear reactions (e is the nuclear energy production rate per unit 



mass), the heating due to turbulence V-F = — V- /jT||_D|| (g> VS 



where |D|| is the 

cddy-diffusivity tensor and the ohmic heating J = \\r]\\ ® A ifj ■ A ifj . 
Its expansion in spherical functions leads to the following equation for the trans- 
port of the temperature fluctuation: 



CB? d* i + $ dIn HAi+ Ep 



dt 



dt 



L(r) 
M(r) 



Mr) + U 4 



Jj_ 

T> 
in 



(3.15) 



T[(r) being given by: 



where 



% = %, v (r) + 7I ;Th (r) + 7I ;N _ G (r) 



(3.16) 



fv (M) _ (e + 
4nGp ( 



+ 



fv.i (M) 

4irGp 
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fv (M) 



t-f cm 
-<Z;Th — 



^drXt (r) 



1(1 + 1) H T 
3r 



and T; ;N _ g = ^ + £ ^> [ Xl ( r ) + (/ e£T - / e + + (/ eCp + / e p - p)Aj] 

with X ; (r) = H T d r ^i - (1 - S + xt)^i - (f + X^)^i- L is the luminosity of the 
star, M its mass and K = x/~pC p the thermal diffusivity. The definition of V, 
Vad, C p , Hp, H T , fe, e, egrav, £t, Xt, £/*, X/ij e m and /°m are given in Mathis 
& Zahn 2004-2005. Terms /-p, /p,; and gi, the fluctuation of the gravity on an 
isobar, are associated to the meridional perturbing force, J--p, namely the sum of 
the centrifugal force, Tc = \Q 2 V (r 2 sin 2 9) , and of the meridional Lorentz force, 
3~c,p- If we project T-p on the vectorial spherical harmonics, its explicit expansion 
in function of f2, £j and xi is derived (cf. Mathis & Zahn 2005). The entropy of 
mixing is taken into account due to the multi-species caracteristic of the stellar 
plasma (cf. Maeder & Zahn 1998), $ being a function of the metal mass fraction 
Z and of ~p. Finally, ohmic heating has been expanded in spherical functions, Ji 
being the associated modal radial functions (their explicit expression in function 
of £z and Xi is given in Mathis & Zahn 2005). 

Terms I, II and III in eq. 13.151 correspond respectively to the multi-species 
characteristic of the stellar plasma, to the advection of entropy by the meridional 
circulation and to the ohmic heating. Next, 7/ : p is the barotropic term of the 
meridional circulation which is generated by the perturbation of the thermal im- 
balance by Tv- ^I;Th corresponds to the thermal diffusion. In fact, if we keep 
its higher-order derivatives, we obtain £s^^d r (Hxdr^i) which is directly related 
with the temperature laplacian. Finally, 7/ ; n-g is the term related to the nuclear 
reactions and to the heating due to the radial adjustements of the star during its 
evolution. 



Finally, the equation for the transport of chemicals: 



d t Ci 



V ■ V c 



= V- (p\\D 



Vc, 



(3.17) 



has to be solved to study the mixing; Ci is the concentration of the i th element 
which is considered. Expanding it on an isobar, it is obtained for the average of 
the concentration of each chemical, cc 



d c t + \d r [r 2 pc l U^] 



1 



At 



3 r [r 2 p(D v + D eS )d r Ci\ , 



(3.18) 



where U^ lS is the velocity associated to the microscopic diffusion processes while 
the strong horizontal turbulence leads to the erosion of the advective transport 
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which becomes a diffusive process (cf. Chaboyer & Zahn 1992) with the following 
effective diffusion coefficient: D c s = S/>o i(i+i)(2i+i)D, • 

Next, taking the definition of the mean molecular weight: 
i = J2i [(1 + %i) /Ai] Ci (Ai and Zi are respectively the number of nucleons and of 
protons of the i th element which is considered), the following advection-diffusion 
equation is obtained for A;: 

dAi din 71, Ui 1(1 + 1) , 

M -A ; - -l-Vp = - V \ DhAi (3.19) 



dt dt H n 



where V„ = f&£. 



Detailed boundary conditions associated to this complete set of partial differ- 
ential equations is given in Mathis & Zahn (2004) and in Mathis & Zahn (2005). 

These equations are all of advection-diffusion type. It has to be underlined 
that advection can not be treated as a diffusion process since it has been shown 
that solving the equation of dynamics begining with a flat rotation profile leads 
to a non-zero differential rotation, the advection being building angular velocity 
gradients (cf. Meynet & Maeder 2000). 

Moreover, note that making hydrostatic and thermal equilibrium hypothesis 
in 13 . 1 II and 13 . 141 returns to the standard stellar evolution equations. 



4 Numerical simulation of secular transport: the hydrodynamical case 
with a 'shellular' rotation 

The numerical simulations presented in this work has been computed with the 
dynamical stellar evolution code STAREVOL and the reader is referred to Siess et 
al. 2000, Palacios et al. 2003 and Palacios et al. 2006 for its detailed description. 
In the hydrodynamical case (B = 0) where we assume that the differential rota- 
tion is shellular, SI (r, 6) = £1 (r), due to the stronger horizontal turbulence which 
enforces the angular velocity to be constant on an isobar, the system of transport 
equations is reduced to: 

(r 2 H) - -^dr (pr 4 nU 2 ) = \d r (p^drTT) (4.1) 
(XV uV T 

and to the I — 2 mode of Eqs. I3.13[ 13.151 and 13.191 (here internal waves are not 
taken into account). It has been now implemented in STAREVOL. Associated 
boundary conditions which are used to compute simulations are detailed in Pala- 
cios et al. (2003). 

The results presented here are issued from the numerical simulation of the evo- 
lution of a 1.5M Q star with a solar metallicity (Z — 0.02) and an initial equatorial 
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rotation velocity VS„i = 100 km.s -1 . The age of the presented model is 7.604 x 10 8 
yr with a central Hydrogen mass fraction X c — 0.57. It is now possible to follow 
for each time-step the internal hydrodynamics of the radiation zone(s) of the star 
which is studied, following simultaneously the differential rotation profile (see Fig. 
Q3-A), the associated meridional circulation pattern (see Fig. [T|-B) and the temper- 
ature and the mean molecular weight fluctuations (see respectively Fig. QJC and 

ED). 




Fig. 1. A: The differential rotation profile. B: Meridional circulation currents. In this 
model, the outer cell is turning counterclockwise allowing the equatorial extraction of 
angular momentum by the wind. C: The T-excesses T^f^ (cos#). It reachs +2.81 x 
10 5 K in the inner region closer to the polar axis and —5.63 x 10 5 K in the inner equatorial 
region and become smaller near the surface. D: The /i-excesses JIA2P2 (cos 9). They only 
occur close to the core and they are positive near the equatorial plane. Dashed lines 
regions correspond to the convective regions. 

Moreover, diagnosis tools have been developed to identify what are the dom- 
inant processes in the angular momentum transport, the meridional circulation 
and the chemicals mixing. 

Considering Fig. [2}A, one can easily identify that here angular momentum trans- 
port is dominated by the advection by the meridional circulation, its flux trans- 
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ported by the shear-induced turbulence being smaller at least by an order of mag- 
nitude except in the region near the center. Then, looking at FigfJJdB, it can be 
easily identified that meridional circulation is mainly driven by the barotropic and 
the thermal diffusion terms. The term due to the nuclear energy production has to 
be taken into account only in the region of the star where nuclear reactions occur, 
here in the center, while the non-stationary term is completely negligible since 
the star which is studied here is a main-sequence star where structural adjuste- 
ments are weak. Finally, concerning the transport coefficients, it can be seen that 
the meridional circulation is the dominating process in the transport of chemicals 
while our fundamental hypothesis concerning the turbulent transport is verified 
(D v «D h ) (cf. Fig.SC). 




Fig. 2. A: Logarithm of the total flux of angular momentum (dots) and of that trans- 
ported by meridional circulation, Fmc (f) = -^ r \pr i QXJ2, (continuous line) and by shear- 

_ 

induced turbulence, Fs (r) = -i- pv v r 4 d r Q, (dashed line). B: Logarithm of the meridional 

circulation (dashed line), the barotropic term (dots), the thermal diffusion term (long 
dashed line), the nuclear energy production and heating due to gravitational adjustemcnts 
term (continuous line) and the non-stationarity term (dashed dotted line) profiles. C: 
Logarithm of the themal diffusivity (dashed dotted line), the horizontal eddy-diffusivity 
(long dashed line), the effective diffusivity associated to meridional circulation (dashed 
line), the vertical eddy-diffusivity (dots) and the total vertical diffusivity, Dt = D v + D e g, 
(continuous line) profiles. Dashed lines regions correspond to the convective regions. 
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Work is now in progress to implement the differential rotation in latitude and 
the transport by the magnetic field and the gravito-inertial waves, those two last 
processes being crucial to understand the internal transport of angular momentum 
in the Sun and the properties of low-masses stars. 

5 Asteroseismic diagnosis 

Since theoretical results which contain some parametrization related to the internal 
dynamical processes are implemented in a stellar evolution code, constraints have 
to be obtained from observations on stellar interiors. Asteroseismology is the 
best way to get such constraints. In fact, rotational transports affects the stellar 
structure during the evolution of the star and thus oscillations frequencies are 
modified. If we recall the classical asymptotic relation for the frequencies of p- 
modes: 

v n ,\ ~ (n + 1/2 + e) Av- ^-tll fo/ 2 , (5.1) 

6 

where Av and 5ij + 2 are respectively the large and the small separation: 
(~ /"^drV 1 , e r„ ,s Ai/ f R dc s dr 

\ JO c s I *~ v n,l JO Q ' ' 

differences are obtained if we consider a non-rotating or a rotating model. More- 
over, constraints could be obtained for each specific transport process by studying 
the effect of its different modelling which could be done on seismic observables. 
For example, if models are computed using the different prescriptions which ex- 
ist for the horizontal eddy- viscosity, Vh, the effect on the small separation once 
the star has evolved is important with an increase associated with that of the 
mean rotation rate of the star and of Vh which leads to a more important trans- 
port and thus to a stronger rotational transport's indirect effect on frequencies. 
Here, results are reported for a 1M© star with different initial rotation velocities 
, Vi = 0, 10, 30, 100 km.s" 1 at t = 8 Gyr (cf. Fig. [3]). The models are built using 
the solar calibration and the different prescriptions for the horizontal turbulence 
established by Zahn (1992), Maeder (2003) and Mathis et al. (2004). Simulations 
are performed using the Geneva stellar evolution code (cf. Meynet & Maeder 
2000) where the microscopic diffusion has been implemented in the same way that 
in Richard et al. (1996), the braking being treated using the formalism of Kawaler 
(1988). The asteroseismic analysis is performed using the Aarhus oscillation code 
(cf. Christensen-Dalsgaard 1997). 

6 Conclusion 

In this work, a coherent description of the dynamical transport processes which 
take place in stellar radiation zones has been undertaken. Each of them and their 
respective effects on angular momentum and chemical transport has been identi- 
fied and modelled in a consistent way. For the first time, a two-dimensional picture 
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Fig. 3. Value of the small difference So,2 in function of the frequency v for different values 
of the initial rotation velocity, Vi for a IMq star. A: Case of the prescription of i/j, of 
Zahn (1992). B: Case of the prescription for vh of Maeder (2003). (From Eggenberger 
Ph. D. Thesis) 



of the internal dynamics of stellar radiation zones is obtained and the first step of 
the numerical implementation of the theoretical results in stellar evolution codes, 
namely the purely hydrodynamical case where it is assumed that the strong hor- 
izontal turbulent transport enforces a shcllular rotation law, has been achieved. 
Moreover, work is in progress to implement the differential rotation in latitude, 
the magnetic field, the gravito-inertial waves and the tides while theoretical work 
is engaged to derive prescriptions for the MHD instabilities, the waves excitation 
and the tides dissipation. On the other hand, we have shown how classical seismic 
observables are affected by the internal dynamical processes that will allow to put 
constraints on our modelling from the comparison with the observations. Thus, 
with the first light of COROT, the preparation of KEPLER, SOHO and forth- 
coming hclioseismic space instrument such as GOLF-NG on DynaMICCS project 
(Turck-Chieze et al. 2005, Turck-Chieze et al. 2006), we are now entering in a 
new exciting period of the story of stellar evolution where we hope to be able to 
obtain a dynamical vision of the Hertzsprung-Russel diagram with the support of 
helio and asteroseismology. 
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